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The general family of kinetic equations, which in plasma kinetics are called the BBGKY equations, 
are obtained rigorously from basic probabilistic considerations in order to exhibit explicitly the con- 
ditions or assumptions under which they obtain. 
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1. Introduction 

The plasma kinetic equations, referred to as the BBGKY hierarchy, are fundamental in the 
analysis of plasma kinetics. In the following, these equations are derived rigorously in order to 
determine explicitly the conditions or assumptions under which they are valid. The heuristic 
derivations do not clearly indicate the assumptions which are important to both theoretical and 
experimental plasma research. Since the functions in the BBGKY equations are probability density 
functions the following treatment is probabilistic. 

2. Preliminaries 

Consider first Euclidean 6-space E e in which a point is represented by the 6-tuple of real 
numbers (or coordinates) {qi, q-i, q-.\, pi, p>>, p-,{\ and n points are represented by n 6-tuples {q s x , 
Q'h Qv Ph Ph Ps) for 5= 1 ... n or more concisely for each 5= 1 ... n by {q- p-} where i= 1,2,3. 
Equivalently these n points may be represented in Euclidean 6/i-space E 6n by one point, or 6n-tuple, 
{qU Pf} where s= 1 • • • ft>i = l»2, 3. 

Consider also a set of 6n elements Q\ . . . Qz n P\ • . . Psn called observables; the entire 
set Qi . . . QznPi • • • P$n being called the System S of observables. We then make the definition: 
Definition: A (static) state P of the System S is an assignment P : Qi . . . QsnPi • • • P*n~* i^l 
. . . q", pj . . . p"} i.e., an assignment of6n numerical values to the observables or an assignment 
of a point in E 6n . The correspondence between the P's and the points {q% p S j} is assumed to be 1:1 
so that we may write P= {q}, p?}. 

The statistics of the System (in this static case) are given by a probability distribution F on 
the Borel Field of E 6n such that the following assumption holds: 

ASSUMPTION 1: F is absolutely continuous with respect to Lebesgue measure and thus there exists 
a measurable function p = p(q?, p?); i.e., a probability density function, such that for any set B in 
the Borel Field one has: 



F(B) 



= P(q' , Pf)(dq dp) n where (dq dp) n = dq} . . . dqj dpj . . . dpj. 
Jb 
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In order to introduce dynamics into this schema, we specify the state concept in the dynamic 
case as follows: 

DEFINITION: A state function P of the System S is a mapping from the direct product S X R # , 
where R # is the real numbers, into E 6n such that for each teR # , P provides an assignment P(t) of a 
point in E 6n to the observahles in S, that is, for each t, P specifies a {static) state P(t) : Qi . . . Q 3n 
Pi . . . Pan - * {qf(t), pf(t)} which by the preceding convention we may write P(t)={q-(t), p?(t)}. 

We shall also make the following assumptions: 
ASSUMPTION 2: For given t and arbitrary {q\, Pj}eE 6n there exists a state function P such that 
P(t)={^(t),pf(t)} = {^,pf}. 

ASSUMPTION 3: The dynamics of the System S are given by a Hamiltoni an function H n = H n ({qj\ p?}) 
and by a one-parameter group {U t }, teR # , of contact transformations, differentiate with respect 
to t, so that ifP is an arbitrary given state function then for teR # , reR # 

U T P(t) = U T {q?(t), p?(t)} = {q?(t + T), p?(t+r)}=P(t + T). 

The statistics of the dynamical System are now given, for each value of the parameter teR *, 
by a probability, distribution F t which for each teR* satisfies assumption 1 above so that we may 
write: F t (B) = p,(<#, pf){dqdp) n . 

JB 

Definition: The System S in the state P(t) will be said to be in the Borel set B at t if P(t) = 
{q s (t), p s (t)}eB. In these terms F t (B) may be referred to as the probability that at teR # the System 
will be in a state in B. 

Since the variation of the states P(t) , or of the state function P, with t is given by assumption 3, 
it is possible to describe the related variation of the F f (B) with t as follows. By assumption 3: 
U T P(t) = P(t-\-r), the states P(t) which are in a Borel set B at t will become states P(t-\-r) in a 
Borel set B 7 at t + r. The set B T is the set of all points {qf , pf} of £ 6 " which satisfy {qf , pf}= {qf(t + r), 
pf(t + r)}= U T {qf(t), pf(t)} for {qf(t), pf(t)} = {q-, pf}eB. Since the U t are contact transformations 
and since the Poincare Invariants under contact transformations include the Lebesgue measure 
fji(B) of the Borel sets B C E 6 \ we have /jl(B t ) = p,(B). 

For given t and Borel set B with state {qf(t), pf(t)} = {qf, pf}eB consider pdqf, pf)/m(B) = 
p(qf(t), pf(t), t)pi(B) and similarly for r and B 7 with {q?{t + r), pf(t + r)} = {qf, p]}eB T consider 
Pt + MU P?)/*(£ T ) =p(q?(t + T),pf(t + T) )fi(B 7 ). 
ASSUMPTION 4: The functions p t (qi\ Pi) = p(q?(t), pf(t), t) are differentiate with respect to q?, p- 

and t and -t- (which exists by assumed differentiability of the contact transformations U t ) is uni- 
formly continuous over each Borel set B. 

Then for t= At, sufficiently small, and arbitrary closed bounded Borel set /3 C B one can write 
by the law of the Mean: 

p(qf(t + At), pf(t + At),t-hAt) / x(/3 At )=p(qf(t + At),p?(t + At),t + At) M (/3) 

= p(qJ(t),p^(t),t) i Lt( i 8)+At^(qf(t + ^At),p^(t + ^At),t + ft,At)/x(/3) 

. . r . . r dp 

and using uniform continuity ol —r one can write 

dt 

F t+At (B At )-F t (B) f dp, s s , (A , . 

, . dFt(B) i i r i i • i 

Letting — — represent the lelt hand expression we make 
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ASSUMPTION 5: — k = Ofor all Borel sets B. From this assumption and continuity of-f- we have 






i=l s =l x M i 



r =-^-+ [p, H n ]. [Liouville's Theorem] 

at 

3. Derivation of BBGKY Equations 

Consider a subspace E a C E 6n where E a is the subspace of 6cr-tuples {qf, pf}, (0=71,72 • • -7o-, 
the 7a being cr distinct integers from 1 to n; for purposes of simplicity in notation take o> = l,2, . . . cr. 

For B^ a Borel subset of E " let C a be the 6ra dimensional "cylinder" on B a i.e., C (T = {{qh P?K 
i = l, 2, 3, s= 1 . . . n; such that for 5 = 1 . . . a {g?, p}*} C B fT }. We shall write C a = B a X E'^ 
where E[ r = {{q^ p?} , 5 = <r+ 1 , . . . n}. 

Since C a is a Borel subset of £ 6/ ' we have as above 

Ft(C a )= p ( (ql p\){dqdp) n 

J(\r 

= \ I Pt(qf, Pf)(dqr x • • • dpi')] (dq> . . . dp?). 
Writing pAql Pi t)= p t (qf, pf)(dqr' . • ■ dpf) 

Ft(C (T )=\ P(T (qlPlt)(dq} . . . dpf). 



th 



The Hamiltonian H n = H n ({q], pi}) will be of the form in the following assumption: 
Assumption 6:. H n = (i] 2 T^") + v ({<Ii}) 

where V({qf}) = f) £ fl k , t + J [^({qf}) + n«q?})l 

k=l l>k k=l L J 



w;^/i k i = 6^i k = ^(q^, q| ) = a symmetric function of q^, q^; for example: kl = e 2 



3 1-1/: 

Kqf-qi) 2 

i = l J 



and 3>w({Q'i}) = potential at kth point (or particle) due to the boundaries (in E 6 ) which are assumed 
to contain the n points (or particles) and y ^+({q s i }) = potential at the kth point (or particle, or elec- 
tron) due to a fixed background electric field (produced by + charges). 



One can then write 
dF l 



f^-JJ,^ <«<'>.*«.***>. 
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-/JJH,(5®-*(5))j(-i*)- 

jB <r Jt cr k = (T+ l W»' \>* °*« m<; 



k 

i_ 

m 



dOnac | d<P k w | d¥* 

a^f a<7f a<?f 



In the above and following we omit the (dg, dp) notations whenever they are obvious. 
Since 



w =0, fc and £ 7^+ 2 



^ agfcj 



f* dqf ^ dqf f-> bq\ 

l>k * l m<k ^ l l=\ ^ l 



the above may be written, using the definition of p a : 









Ij,. fc 5 +1 (4)l54' 






and one can abbreviate the notation for these five integrals as: 
=7(1) -7(2) -7(3) +7(4) -7(5). 

Now in the above notation for E a = {{q k \ Pi'}* k ~ 1 • • • °"} letting E {1) = {{q k , p k }; k = 1} we 
set E a ' l = E (r + E in , that is, the minimum subspace containing the subspaces E a , E il) . Also, as 
above, let E' (T = the complementary subspace of Z? 47 in E 6n consisting of all 6(n— o-)-tuples {q k , p k } 
where k = cr+ 1, crH-2, . . . ti and i= 1, 2, 3 and let £^ , = the complementary subspace of E a < l 
in E 6n consisting of all 6(n— a— l)-tuples {q k , p k ) where ^/, cr + 1 ^ A: ^ n and i = 1, 2, 3. Also 
similar to the definition of p a above we shall let: 



P<r,i(QU Ph ')= p(Qh Ph t)(dqdp)' a ^ 



where (dq dp)' a ^ represents the usual "volume" element in E^ ,. 
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Furthermore we now make the following: 
Assumption 7: p(q s , p s , t) is a function symmetric with respect to the indices s. 

It is clear then that the functions p (T+ \(q > j, pf, t) =p<r, o-+i(<7?, PU l ) an ^ Pa, i(#f, pf, differ 
only in that the coordinates qf +1 , pf +] occurring in p a +\ are replaced by q\, p\ in p^, /. Thus 

J B o- k=l *i ^k* 7 L(r M i 

-/„| 1 ^ 1 |/ E0) l|L./K.p;..)w< I dP);.,w<dp:) 



since 






Suppose now that the states of the System are restricted by: 
ASSUMPTION 8: There exists a closed bounded subset 2 3n o/E 3n such that the probability density 
function p(q-, p-, t) is identically for all (q?, p?, t) such that {q\} s=l . . . n,i=l,2,3wa 
point on the boundary of 2 3n or outside £ 3n . This is equivalent to saying that for any t and any 
Borel set B with B C\ E 3n in the closure of the complement o/£ 3n , the probability, F t (B), that the 
System is in a state in B is 0. 

We shall also take the boundary of 2 3n to be the boundary pertaining to the potential functions 
Og, of assumption 6. 

One can then write: 

Too QQ 

and I — *j dqf = by virtue of assumption 8 hence 7(4) =0. 

If also one makes: 
Assumption 9: For any giVe/i k—1 . . . n p(q-, p-, t)— » as p*^±oc 
then 1(5) =0 by the same type of argument. 

Finally from assumption 5 — -3 — — =0 for any C a hence for any B a = 1(1) —1(2) — 1(3) 
and since the integrands are assumed continuous we get the BBGKY hierarchy equations: 

2£z+ v (*£z\ ^-V ^ (x> ddkl 1 ^ 1 d * k A 
dt "*" 2, Wv W 2/, VM7 V,- 2 ? atf art dtf 7 
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If one takes the Hamiltonian H a = H (T {{q s i , p?})s = 1 ... cr to be 
then these equations become 



?+wj-I*— >/„ s -fr ^? <*•■*? 



The conditions under which these BBGKY equations are valid are thus explicitly given by the 
assumptions 1 through 9. 

4. More General Case 

The preceding discussion is conditioned by its restriction to Systems of 6n observables 
Qi, Pf i— 1 . . . 3rc, that is, to exactly n points (particles) in E 6 . Therefore the functions 
F r {B) — F(B, t) are, in this sense, conditional probability distributions. A somewhat more general 
situation may be treated as outlined briefly in the following. 

Let E™ be the linear space of all infinite double sequences: 

{q\, qh qh qh qh • • • Pi> Ph Ph Pi> Ph • • .} = {??, pf} * = 1, 2, 3 s=l, 2, . . .. 

If E G = {qt, pi} & = 1, 2, 3 then E w may also be regarded as the infinite direct product: 

E" = E 6 xE 6 xE 6 x. . .. 

Similarly writing as above E a = {{<$, pf}; i = l, 2, 3; 5 = 1, 2 . . . cr} and B a for a Borel subset of E a 
then for m ^ cr we denote the cylinders Cg? in E m and C% in E™ by: C™ = B a XE™' where 

E?'={{qi,Pth s = <r+l 9 . . . m} and C%=B<rXE%' where E%'={{ql pf}; s = cr+l, . . .}. 

If the System may vary stochastically with respect to t we may let P m (t) be the probability 
that, for teR , there are m points (particles) in E 6 or equivalently that the System consists of 6m 
observables Qi, Pu Also let F™{C™, t) be the conditional probability that, if there are 6m observables 
in the System, then this System will be in a state P={qf, pf} such that the 6cr-tuple 
{#!' Pf} s=l ... a will be in B a at t. Then the probability that the System will be in a 

state P= {#, pf} which is in B (T at * is given by: F%{C%, t)= ^ P m (t)F™(C™, t). 

If for each m the distribution F™(C™, t) is absolutely continuous with respect to Lebesgue 
measure in E m then there exists as above a measurable function pJJ! such that 

F*(C* = f p»'(<?f, pf, (dq dp) m . 
Jc a 

Furthermore F^(C£, t) will also be absolutely continuous with respect to the product measure E u 
and F%{C a , t) = I p%(q\ ph t) {dq dp)^ where p% is a function measurable with respect to the 
product measure in E™. 

n 

Since F%(C%, 0=2 p m{t)F%(C™, t) the sequence j P m (t)p^(ql pf, t), as rc^oo, C on- 
verges in measure to p£(gf, pf, and under proper assumptions on the uniform continuity of the 
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derivatives of the terms of the sequence one may relate the hierarchy equations for the p£ to the 
above derived corresponding relations for the individual pJJ? occurring in the terms of the sequence. 

5. Summary 

The preceding derivation of the BBGKY hierarchy equations and the specific assumptions on 
which it is based are not restricted to plasma kinetics. As indicated in the introduction the treat- 
ment is essentially probabilistic and nearly all the assumptions are probabilistic in nature. The 
excepted nonprobabilistic assumptions refer to the existence of a Hamiltonian and an associated 
contact transformation group, both of which are more general than the requirements of plasma 
kinetics. 
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